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The variety of 7-ray burst phenomenology could be largely attributable to differences in the 
opening angle of an isotropic outflow or to a standard type of event viewed from different 
orientations. Motivated by this currently popular idea, we study the effects of varying the 
energy per unit solid angle in the unsteady expelled outflow by an increase either in the 
bulk Lorentz factor or in the baryon loading. We apply these models to interpret the ob- 
served correlations between variability, luminosity and spectral peak energy and find that 
while the latter scenario fails to provide a good description, bursts produced by collisions 
between similar mass shells but with increasingly large Lorentz factors are both more vari- 
able and have larger peak spectral energies. We present detailed internal shock calculations 
confirming this interpretation and discuss the roles various timescales, radii and the optical 
thickness of the wind play in determining this wide range of behaviors. Finally, we dis- 
cuss the variety of scenarios in which large variations of the source expansion velocity are 
naturally expected. 
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1 Introduction 



The discovery of afterglows in recent years (Costa et al. 1997; van Paradijs et al. 
1997; Frail et al. 1997) has moved the study of gamma-ray bursts (GRB) to a new 
plane. It not only has extended observations to longer timescales and other wave- 
length bands, enabling the measurement of redshift distances, but has also provided 
a direct determination of the source size and confirmation of its relativistic expan- 
sion (see Piran 1999 and Meszaros 2001 for recent reviews). Detailed predictions 
of the afterglow properties, made in advance of the observations (Meszaros & Rees 
1997; Vietri 1997), agreed well with subsequent detections at these wavelengths, 
followed over timescales of months. While progress has been made in understand- 
ing how the GRB and afterglow radiation arises in terms of an unsteady relativis- 
tic outflow followed by the development of a blast wave moving into the external 
medium, interest continues to grow as new observations provided new challenges 
of interpretation. 

Many GRBs are now found to be associated with star forming regions (Kulkarni 
et al. 1998; Fructher et al. 1999; Berger, Kulkarni & Frail 2001), the remnants of 
a massive stellar progenitor system (Piro et al. 2000; Amati et al. 2000) and possi- 
bly supernovae (Galama et al. 1998; Bloom et al. 1999; Reichart 1999; Bjornsson 
et al. 2001; Lazzati et al. 2001b). These observations give support to the idea that 
the most common GRBs could be linked to the collapse of massive stars (Woosley 
1993; Paczyhski 1998; MacFadyen & Woosley 1999). This speculation is made 
more intriguing by the recent report that GRBs could be highly collimated, so that 
there may be many more events Doppler boosted in directions other than ours, 
bringing the GRB rate to within a factor of 100 or so of the supernova rate and their 
derived explosive bulk energies down to ps 10 525 erg (Frail et al. 2001, Panaitescu 
& Kumar 2001; Piran et al. 2001; assuming a conversion efficiency of 1%). The ap- 
parent dispersion in isotropic energies is then caused by a distribution of jet opening 
angles - bursts with higher energy per solid angle are more collimated (or maybe 
viewed closer to the jet axis; see Rossi, Lazzati & Rees 2002; Salmonson & Galama 
2002; Zhang & Meszaros 2002). New evidence and new puzzles were added when 
the total isotropic luminosity was found to be also related with a variety of tempo- 
ral and spectral behaviors exhibited by 7-ray light curves - namely the degree of 
variability or "spikiness" (Fenimore & Ramirez-Ruiz 2002a), the differential time 
lags for the arrival of burst pulses at different energies (Norris, Marani & Bon- 
nell 2000, Norris 2002), and the rest frame GRB peak energy (Lloyd-Ronning & 
Ramirez-Ruiz 2002). 

These correlations are still tentative, but if confirmed they could be used to relate 
the properties of material moving almost directly towards us (e.g. the variable ac- 
tivity of the 7-ray light curve) with the ejection in directions away from our line 
of sight (e.g. degree of collimation; see Salmonson 2000; Ioka & Nakamura 2001; 
Kobayashi et al. 2002; Plaga 2001). This could turn out to be very important, since 
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at present we can only infer the energy per solid angle; as yet the constraints on 
the angle-integrated 7-ray energy are not strong. The gamma-rays we receive come 
only from material whose motion is directed within 1 /Y of our line of sight. They 
therefore provide no information about the ejecta in other directions. At observer 
times of more than a week, the beaming and aberration effects are less extreme so 
we observe afterglow emission from a wider range of angles. At these later times, 
however, all memory of the initial time-structure of the outflow responsible for the 
7-ray emission would be lost. 

Information regarding how much energy and momentum has been injected, its dis- 
tribution in angle, the mass fractions in shells with different Lorentz factors and its 
efficiency in radiated gamma-rays are at the forefront of attention. We address some 
aspects of these issues here by studying the physical conditions under which sim- 
ulated GRBs arising from internal shocks in relativistic winds can account for this 
burst phenomenology - degree of collimation, variable 7-ray activity, the rest frame 
peak energy and the total burst luminosity. To that effect, we consider anisotropic 
relativistic outflows where changes in the energy per solid angle are caused either 
by an increase in the bulk Lorentz factor or in the baryon loading. We apply these 
models to interpret the large variety of behaviors exhibited by both the afterglow 
and the prompt emission, and discuss their possible use for predicting the internal 
structure of the collimated outflow. Some possible scenarios in which large vari- 
ations of the energy per solid angle are naturally expected are highlighted, along 
with the types of observation that would discriminate among the various models. 
We assume H = 65 kms" 1 Mpc" 1 , £l matt er = 0.3, and Q A = 0.7. 



2 Anisotropic or Highly Beamed Outflows? 

The high Lorentz factors and energies seen in GRBs are consistent with the catas- 
trophic formation of a stellar black hole of a few M , with 1 % going to a 
relativistic outflow. This could be the extreme example of the asymmetric explo- 
sion produced by supernova (Khokhlov et al. 1999), in which instead of halting at 
the neutron star stage, the collapse continues to the black-hole stage, producing an 
even faster jet in the process (MacFadyen & Woosley 1999). GRBs arising from a 
very small fraction of stars that undergo this type of catastrophic energy release are 
likely to produce collimated outflows. 

Even if the outflow is not highly collimated, some beaming is expected because en- 
ergy would channeled preferentially along the rotation axis. Also, one would expect 
baryon contamination to be lowest near the axis, because angular momentum flings 
material away from the axis and material with low-angular momentum falls into 
the black hole. The dynamics, however, are complex. While numerical simulations 
of collapse scenarios can address the fate of the bulk of the matter (MacFadyen & 
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Woosley 1999; MacFadyen, Woosley & Heger 2001; Aloy et al. 2000; Zhang et al. 
2002, in preparation), higher resolution simulations of the outer layers of the stellar 
mantle seem to be required since even a very small amount of baryons polluting 
the outflow could severely limit the attainable Lorentz factor. It is quite possible, 
for instance, that the stellar pressure will tend to collimate the fireball into a jet. 
A broad spread of Lorentz factors is thus expected - close to the rotation axis T 
may be high. At larger angles away from the jet axis, there may be an increasing 
degree of entrainment, with a corresponding decrease in T. So in this model, the 
measured flux is more intense when observed closer to the jet axis (with an angle 
9 V pa from the jet axis). A large spread in the inferred luminosities thus origi- 
nates from the differences in viewing angles (see Figure la). This interpretation 
may be appealing because a large spread of both Lorentz factors and luminosities 
could be attributable to a standard type of event viewed from different orientations 
(Rees 1999). On the other hand, the large variations in the fluences of GRB could 
also originate from differences in the opening angle 6j of an isotropic outflow (see 
Figure lb; Frail et al. 2001; Panaitescu & Kumar 2001; Piran et al. 2001). The dy- 
namics of the blast- wave in the off-axis anisotropic jet are almost indistinguishable 
from the uniform jet model, provided that the energy per unit solid angle varies as 
6~r 2 (Kobayashi et al. 2002; Rossi et al. 2002; Salmonson & Galama 2002; Zhang 
& Meszaros 2002). 

The correlations between variability, the luminosity per unit solid angle and the 
characteristic photon energy in the cosmological rest frame (Lloyd-Ronning & 
Ramirez-Ruiz 2002) impose severe restrictions on any emission model. The ques- 
tion at hand is whether this burst phenomenology, combined with the intrinsic vari- 
ations in energy per unit solid angle discussed above, can be understood in terms of 
internal shocks and if so, under what conditions. To this end, the results presented 
in §3 are independent of whether the variations in energy per unit solid angle are 
caused either by the viewing angle of a standard event (provided that dT/d6 v is 
small for dO v < T" 1 ) or by collimated jets, in which the typical Lorentz factor (or 
energy per unit solid angle) is higher when the jet is highly collimated. 



2. 1 Highly collimated outflows 

As discussed above, the gamma-rays we receive come only from material whose 
motion is directed within an angle T" 1 of our line of sight. At observer times of more 
than a week, the shocked outflow starts to slow down. As V drops after the decel- 
eration shock, the causal angle includes an increasing amount of the solid angle 
toward the jet as well as towards the equator, so we observe emission from a wider 
range of angles. The afterglow is thus a direct probe of the geometry of the ejecta - 
if the outflow is beamed, we expect a downturn in the light curve as the edge of the 
jet becomes visible (Rhoads 1997). Collimation factors of Vlj/Aix < 0.01 have been 
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derived from the observation of such steepeningsQ (Kulkarni et al. 1999; Castro- 
Tirado et al. 1999; Harrison et al. 1999; Frail et al. 2001; Panaitescu & Kumar 
2001). Although several methods of producing and collimating jets have been pro- 
posed, these highly collimated outflows of plasma with velocities close to the speed 
of light may be a unique feature of rapidly rotating, gravitationally confined plas- 
mas that are threaded by a strong magnetic field (e.g. Uchida & Shibata 1985). An 
important result from present numerical simulations of Poynting flux-dominated 
outflows, which was not predicted by semi-analytical steady state studies, is that 
the jets arising from these outflows are likely to be highly collimated (see Meier, 
Koide & Uchida 2001 for a recent review). Because of the strong pinch that devel- 
ops, a narrow jet that delivers its thrust in a narrow solid angle may be a common 
ingredient of strong rotating magnetic fields, not only for accretion disks but also 
perhaps for many jet-producing objects. 

The range of parameter space where high jet speeds are expected (i.e. jets propelled 
at speeds much larger than the local escape velocity) is determined by the condition 
^mhd > Lc- The MHD luminosity 



depends only on the poloidal magnetic field B p protruding from the jet-production 
region at radius R and on the angular velocity of the field Q (Blandford & Zna- 
jek 1977), while the critical luminosity L c is defined to be that given by the ratio 
between the energy needed for the plasma to reach the escape velocity and the 
free-fall time, 



where p is the density of the plasma. The condition L MH d > L c is a necessary 
requirement to explain the T > 100 outflow seen in GRBs, even from the vicin- 
ity of black holes. This criterion is similar to that of LeBlanc and Wilson 1970, 
in which the magnetic field grows to a dynamically important value in less than a 
dynamical time, which seems to be a condition of forming MHD jets in collapsing 
supernova cores (Meier et al. 1976; Wheeler et al. 1999). Although carrying more 
power, these highly collimated outflows will be much less efficient in imparting 

3 The above argument assumes that the breaks observed in many GRB aftergiow light- 
curves are due to a geometrical beam effect and not to either a transition to non-relativistic 
expansion (Huang, Dai & Lu 2000) or an environmental effect such as a sharp density gra- 
dient (Chevalier & Li 2000; Ramirez-Ruiz et al. 2001). An important feature produce by 
this jet collimation is the achromaticity of the afterglow break, which will clearly distin- 
guishes it from the steepening that may be produced by the passage of a spectral break 
through the observing band. 
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energy and momentum to the outer stellar layers and may not explode (Khokhlov 
et al. 1999). These failed explosions may be similar to the failed supernova model 
for GRBs (MacFadyen & Woosley 1999), continuing to accrete much of the stellar 
mantle and thus collapsing into a black hole. 

While the very narrow 2" and fast jets T > 100 that are speculated to exist in 7- 
ray bursts (Panaitescu & Kumar 2001) may be a general feature of strong rotating 
magnetic fields, some observations may be difficult to reconcile with this interpre- 
tation. For example, several afterglows have shown evidence of large amount of 
X-ray line emitting material, possibly arising from ionized iron (Piro et al. 1999; 
Yoshida et al. 1999; Piro et al. 2000; Antonelli et al. 2000). The large Fe Ka equiv- 
alent widths inferred from the X-ray observations favor models in which the line is 
produced when the primary X-ray emission from the source strikes Thomson-thick 
material and Compton scatters into our line of sight (Lazzati, Campana & Ghisellini 
1999; Bottcher & Fryer 2000; Vietri et al. 2000; Rees & Meszaros 2000; Ballan- 
tyne & Ramirez-Ruiz 2001, Lazzati et al. 2001). These conditions impose strong 
constraints on the location and geometry of the optically thick reflecting material, 
as well as on the structure of the outflow producing the primary X-rays - illumi- 
nation of the stellar remnant may be complicated if the GRB emission is highly 
beamed. This problem may be overcome in particular source geometries for which 
lower luminosities and softer spectra are expected at the edges of the relativistic 
outflow (Ballantyne & Ramirez-Ruiz 2001). More importantly, highly collimated 
outflows bring the GRB rate to within a factor of 100 or so of the supernova rate 
(Frail et al. 2001). Those stars that produce GRBs are likely to have core masses 
> 6M Q at the time they collapse (Woosley 1993). Stars satisfying this criterion are 
pa 30-50 times rarer than those producing SN Ib/c (Izzard, Ramirez-Ruiz & Tout 
2002, in preparation). Because burst formation is also likely to be favored by rapid 
rotation and low metallicity (MacFadyen & Woosley 1999), highly collimated out- 
flows could face some difficulties in explaining large event rates (i.e. 1/100 of 
the supernova rate). On the other hand, one appealing aspect of massive star pro- 
genitors is that the large variety of stellar parameters could explain the diversity of 
jet opening angles. If however the GRB 980425/1998bw association is real, then 
we may have a new class of GRB with lower energy £ 7 « 10 48 (f2 ; /47r) erg, which 
is only rarely observable even though its comoving density could be substantial. 



2.2 Anisotropic relativistic outflows 

Broader jets can also be produced by a magneto rotational mechanism, however, 
some concern has been raised that these outflows may not be easy accelerated to 
T > 10 (Meier et al. 2001 and references therein). Detailed predictions are nonethe- 
less still in need of three dimensional, high-resolution core-collapse calculations. 
Other mechanisms of jet production such as neutrino radiation in the context of 
the collapsar model (MacFadyen & Woosley 1999) or intense radiation of a newly 
born pulsar (Blackman & Yi 1998) could be less constraining in their beaming 
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requirements. 



If the large variety of behaviors exhibited from burst to burst is produced by viewing 
an anisotropic universal beam configuration from different directions, the opening 
angle of the jet 9j > 9 V should still affect the light curve at the time when the edge 
of the jet becomes visible (this break would be additional to the one predicted 
from the viewing angle dependence; Rossi et al. 2002, Zhang & Meszaros 2002). 
If the critical Lorentz factor is less than 3 or so (the opening angle exceeds 20°) 
such a transition might be masked by the transition from ultra relativistic to mildly 
relativistic flow, so quite generally it would be difficult to limit the opening angle of 
this anisotropic outflow if it exceeds 20°. Although, the jet sideways expansion in 
this model is likely to be very important and more detailed calculations are required 
in order to understand the exact shape of the afterglow light curve, it is nonetheless 
important to mention that some afterglows are unbroken power laws for over 100 
days (e.g. GRB 970228), implying that the opening angle of the late-time afterglow 
at long wavelengths is probably > 20°. 

Also, under this interpretation, the 7-rays are likely to be more narrowly beamed 
than the optical afterglow and therefore there should be many "orphan" afterglows. 
The transient sky at faint magnitudes is poorly understood, but there are two ma- 
jor searches to find supernova down to R=23 (Garnavich et al. 1998; Perlmutter 
et al. 1998). These searches are sensitive to afterglows of the brightness levels ob- 
served to date and have covered a few tens of square degree years of exposure. 
It then follows that the afterglow rate is not more than a few times 0.1/sq deg/yr. 
Considering that the magnitude limit of these searches allows optical counterparts 
brighter than 1 ph cm" 2 s" 1 to be detected, the ratio of orphan afterglows to GRB is 
unlikely to exceed ^10 (Meszaros, Rees & Wijers 1999). This condition is not in 
disagreement with the few number of orphan afterglows that are predicted from an 
anisotropic universal beam configuration (Rossi et al. 2002). Studies of the transient 
sky at magnitudes down to R 25 should be able to disprove or consolidate this 
anisotropic model. The location of GRBs within their parental galaxy could help 
settle this question, since recently it was suggested that bursts located closer to the 
center of their parent galaxies have smaller isotropic equivalent energies (Ramirez- 
Ruiz, Lazzati & Blain 2002). If confirmed in further host observations, this corre- 
lation will strongly complicate this universal beam interpretation. A direct test of 
this model remains, however, the time dependent measurements of polarization. 



3 Internal Shocks: A Phenomenological Study 
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3.1 Model Outline 



We simulate GRB light curves by adding pulses radiated in a series of internal 
shocks that occur in a transient, unsteady relativistic wind (Rees & Meszaros 1994). 
Several authors have modeled this process by randomly selecting the initial condi- 
tions at the central site (Kobayashi et al. 1997; Daigne & Mochkovitch 1998; Pilla 
& Loeb 1998; Panaitescu, Spada, Meszaros 1999, Spada, Panaitescu & Meszaros 
2000; Ramirez-Ruiz, Merloni & Rees 2001). The model used here is similar to that 
described by Ramirez-Ruiz et al. (2001) but differs in the following aspects: (i) Our 
treatment of the radiation emission takes into account synchrotron emission and 
Compton scattering; (ii) The effect of photon diffusion through the optically thick 
wind is included, which may be very important in obscuring pulses that occur at 
small radii; (iii) We use the shock jump equations to determine the physical condi- 
tions in the shocked fluid, and (iv) The duration of the pulse width At is calculated 
by adding in quadrature the angular spreading time AT a i?/(2r 2 c), the radiative 
cooling time AT r ss t 7 /T (where t 1 is the lab-frame radiative timescale), the shell 
crossing time ATa ~ t^/T 2 (where t& ~ A/|v sh -v | is the lab-frame shock's cross- 
ing time, v s h is the shell preshock flow velocity and A is the shell thickness) and 
the diffusion time AT d through the optically thick wind as done by Panaitescu et al. 
1999 and Spada et al. 2000. 

The wind is discretized into a sequence ofN = T dm / 5t t shells with a range of initial 
thickness A,-, where rd ur is the duration time of the wind ejection from the central 
source and 5u -C is the average interval between consecutive ejections. The 
time interval 8tj between two consecutive ejections i and i + 1 is assumed to be 
proportional to the i+ 1th shell energy. This implies that longer periods of quies- 
cence, during which the engine accumulates fuel, are followed by more energetic 
shells (this is motivated by the observations of Ramirez-Ruiz & Merloni 2001 and 
Nakar & Piran 2002). The different values of 5t { are selected from the log normal 
distribution of pulse intervals found by Norris et al. (1996). 

We calculate the radii where shells collide and determine the emission features for 
each pulse. If some inner shell moves faster than an outer one (r, > Tj), it will 
overtake the slow one at a radius RiiUj) = RjiUj) = R c oc Tj. For each collision be- 
tween two shells there is a reverse and a forward shock. The shock jump equations 
determine the physical parameters of the shocked fluids: the velocities of the shock 
fronts, the internal energy in the shocked fluid frame u', and both the thickness A,- 7 - 
and Lorentz factor r,y of the merged shell at the end of the collision. We assume 
that in between two consecutive collisions the thickness of the shell increases pro- 
portionally to the fractional increase of its radius ^ oc ^ (e.g. Spada et al. 2000). 
The ejection parameters determine the dynamics of the wind and the dynamical ef- 
ficiency, €ij. This efficiency reflects the differences between the Lorentz factors of a 
pair of colliding shells. The efficiency for an individual collision can be calculated 
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from the initial and final bulk energies, 

e-1- (3) 



where 



r 2 . = rr — LJ. a\ 

J MiTj+MjFi 1 W 



and the resulting mass is M ( - 7 - = Mi+Mj. The first collisions remove the initial 
random differences between the Lorentz factors of successive shells. If the mean 
Lorentz factor f remains steady for the entire burst duration, then the efficiency 
steadily decreases during the wind expansion. If f is modulated on a timescale 
much smaller than the overall duration of the wind, dynamically efficient collisions 
at large radii are still possible. 

A fraction ( < 1 of the shock- accelerated electrons is assumed to have a power 
law distribution -p in the electron Lorentz factor 7 e , starting from a low random 
Lorentz factor 7„, f ^fCH/- 1), where the energy stored in electrons is a fraction 
e e , T- , is the Lorentz factor for internal motions in the shocked frame and we assume 
p = 2.5. The magnetic field is assumed to be turbulent and parameterized through 
the fraction e B of the internal energy it contains: B 2 xs %ne B u' (primed quantities 
are measured in the comoving frame). The synchrotron spectrum of each pulse is 
approximated as three power-law segments, with the high energy slope(s) depend- 
ing on p and the relative values of the peak and cooling frequencies. The typical 
energy of synchrotron photons as well as the synchrotron cooling time depend on 
the Lorentz factor 7 e of the accelerated electrons and on the strength of the mag- 
netic field. The two energy release parameters ( and e B alter the peak energies of 
the synchrotron and inverse Compton spectra (they are higher for lower Q and the 
Comptonization parameter (which is larger for smaller e B ). 

The shock- accelerated electrons radiate, and the emitted photons can be upscattered 
on the hot electrons (7, >> 1) or downscattered by the cold ones (p( e « 1). The opti- 
cal depth to upscattering, far from the Klein-Nishina regime, is Ti c aj h (n' e mm(ct' yl A') 
(Spada et al. 2000), where n' e is the electron density and t' = t'/(l +y) is the ra- 
diative timescale with t' the synchrotron cooling time and y the Comptonization 
parameter (y = 7^ri c for r- lc < 1). A fraction min(l,Ti C ) of the synchrotron pho- 
tons is inverse Compton scattered «i C =max(l, t£) times unless the Klein-Nishina 
regime is reached (Panaitescu et al. 1999). The energy of the upscaterred photons 
is hvi C = min[7 m m,;C 2 ,/if syn (47^/3)" ic ]. The optical thickness r c for the cold elec- 
trons within the emitting shell is calculated by taking into account the cold elec- 
trons within the hot fluid - those that were accelerated but have cooled radiatively 
while the shocked crossed the shell- and those within the yet unshocked part of the 
shell. When r c > 1, photons are downscattered by the cold electrons before they 
escape the emitting shell, leading to a decrease in photon energy and increase in 
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photon duration. We approximate the increase in pulse duration due to the diffu- 
sion through optically thick shells by the time Ar d it takes for a photon to diffuse 
through them (Spada et al. 2000). The addition of all pulses gives both the burst 
spectra and the gamma-ray light curve. The latter is binned on a timescale of 64 ms 
and is used to computed the burst variability, that is, the average mean-square of 
the count variations relative to a smoothed time profile (calculated using a boxcar 
function with a timescale equal to 30% of the duration that contains 90% of the to- 
tal counts; see Fenimore & Ramirez-Ruiz 2002a). A simulated profile depends on 
the specific values of the relevant model parameters describing the wind ejection 
(N 50, St h T max , and r min ) and the radiative efficiency of the pulses (e e , ( and e B ). 

3.2 On the internal structure of the relativistic wind 

We study the effects of varying the the energy per unit solid angle of the expelled 
outflow by an increase either in the bulk Lorentz factor or in the baryon loading. 
One simple possibility is that the central engine ejects consecutive shells in which 
the Lorentz factors T, are randomly selected between r min and r max . The shell mass 
M is then assumed constant and determined by the wind isotropic equivalent energy 
by requiring that Mc 2 = r,-. In this scenario, narrow jets (or jets viewed 

close to the jet axis) should have higher Lorentz factors. Assuming r max oc 6~j 2 , 
then E 4v oc T max (Kobayashi et al. 2002). Nonetheless, it could be also possible 
that changes in the energy per unit solid angle are mainly caused by variations 
in the baryon loading of the wind. In this case, the average bulk Lorentz factor 
T = r min /2 + T max /2 remains unchanged throughout the entire wind and equal to 
T = E 4w /J2f = iMiC 2 with Mi randomly drawn from the interval [M min ,M max ] (the 
shell Lorentz factors T, are random between T-Ar and T+Ar). Assuming a corre- 
lation between M max and 6j, v , then E 4n oc M max . These two extreme scenarios will 
be referred to in the following as impulsive and lethargic outflows, respectively. 

3.3 Impulsive outflows 

To simulate the effects of an impulsive outflow, we calculate the temporal and spec- 
tral profiles of 10 2 realizations for a given r max , which we used to evaluate the 
peak of the spectrum in vF v and the mean variability. The wind (isotropic equiva- 
lent) energy, in all cases, is drawn from a log normal distribution with an average 
value E 4w = 10 54 (r max /10 3 ) erg and a dispersion a E = 10 1 / 3 erg. Numerous colli- 
sions happen during the evolution of the wind. Each collision produces a pulse, 
whose strength strongly depends on T max , since the wind luminosity increases with 
r 

L max- 

The optical thickness r c is determined mainly by the wind luminosity, L w = E 4n /T dm , 
and the range of Lorentz factors in the wind. For small collision radii R oc T 2 , the 
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pulse duration is dominated by the time it takes the photon to propagate through 
all the shells of optical thickness above unity. Thus, photons are downscaterred 
by the cold electrons before they escape the wind, leading to an increase of the 
pulse width and a decrease in the photon energy - photon diffusion through the 
wind widens pulses that would otherwise would appear narrow since they occur at 
small radii (Figure 2). This would be the case when the typical Lorentz factor T, 
is relatively small. The smoothing effect is expected to be stronger when the the 
density of scattering is amplified by new e ± pairs are formed in the originally thick 
scattering medium (Thompson & Madau 2000; Guetta, Spada & Waxman 2001; 
Meszaros, Ramirez-Ruiz & Rees 2001; Beloborodov 2002). For larger radii, the 
pulse duration is dominated by both shell crossing and radiative timescales (Figure 
2), which decrease as T,- increases ( ATa increases as a result of the shell widening, 
while AT r is larger for later collisions because B and 7 m are lower). The latter 
timescale is found to be comparable to the angular spreading one during the whole 
wind expansion, for the assumed linear shell-broadening between collisions. This 
implies that for large Lorentz factors, internal collisions will produce pulses with 
small A7J) and high E 4n - giving rise to a very spiky and luminous light curve. 

By running models with 10 2 different values of r max ranging from 100 to 3 x 10 3 
(assuming r max /r min =20), we have found that more variable profiles (> T max ) have 
higher luminosities and can reasonably fit the observational data (see Figure 3, 
where the shaded region depict the la region). The variability reaches an asymp- 
totic value around V ~ 0. 1 because the light curve is binned on a timescale of 64 
ms to compare with BATSE temporal resolution. On the other hand, for very low 
radii (small T max ), the first collisions remove the initial random differences between 
the Lorentz factors of successive shells before reaching the radius of transparency. 
As a result, the efficiency steadily decreases during the wind expansion giving rise 
to very under luminous burst in which the variability measurement saturates due to 
count statistics. In this particular example, the values of the ejecta Lorentz factors 
are such (i.e. r max /r min is set to be constant) that the average efficiency of indi- 
vidual collisions ejj is independent of the intrinsic luminosity of the outflow (or 
r ma x)- Within our model more efficient collisions can alternatively be produced by 
increasing r max , while leaving r m ; n unchanged. The difference in the shell Lorentz 
factors (i.e. r,-y) in this case will be much larger, giving rise to more luminous pulses 
than in the previous case. A smaller range of T max values is thus required to pro- 
duce similar changes in the burst isotropic luminosity (see the banded region in 
Figure 3). The numerical results shown in Figure 3 assume that the fraction e B is 
distributed uniformly in logarithmic space between 10" 2 and 10" 1 , rism = 1cm" 3 (in 
this case most collisions take place before the deceleration radius of the wind; see 
Fenimore & Ramirez-Ruiz 2002b), e e = 0.25 and that ( m 0.1. 

The process by which the dissipated energy is finally radiated depends on the en- 
ergy distribution of protons and electrons in the shocked material and on the values 
of the comoving density and magnetic field. The internal shocks heat the expanding 
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ejecta, amplifying the preexisting magnetic field or generating a turbulent one, and 
accelerate electrons, leading to synchrotron emission and inverse Compton scat- 
tering. By varying the injection fraction ( one can study the relative intensity of 
the synchrotron and inverse Compton components and determine those values that 
maximize the received flux in the BATSE range. For ( > 10~ 2 the inverse Compton 
emission occurs in the Thomson regime and carries most of the burst energy in the 
BATSE window, while for ( < 10" 3 synchrotron dominates over inverse Compton 
scattering as the latter takes place in the Klein-Nishina regime (Panaitescu et al. 
1999). Figure 4 shows the evolution of the synchrotron (which for (=1 lies mainly 
below the BATSE window) and inverse Compton peaks before and after downscat- 
tering for both thick and a thin winds. For smaller collision radii (i.e. small r max ), 
the emission takes place where the wind is optically thick, leading to a decrease in 
the photon energy - diffusion through the wind decreases the photon energy that 
would otherwise would appear large since they are produced at small radii where 
B is larger. For the set of parameters considered here, the Thompson limit is usu- 
ally a good approximation to treat the downscattering of the synchrotron photons. 
The inverse Compton emission, however, peaks at large comoving frame energies, 
and the general cross section has to be considered. For these energetic photons, the 
cross section strongly depends on the photon energy and varies after each photon 
interaction. Figure 5 shows the shifting of the burst emission towards higher en- 
ergies with increasing variability, due to the increase of r max from 100 to 3 x 10 3 
and assuming r max /r m i n =20 (corresponding to the internal shock parameters of the 
simulations shown in Figure 3). This increase in energy stops once the average bulk 
Lorentz factor is large enough that collisions take place where the wind is optically 
thin (see the shaded region in Figure 5). The typical radiation energy of each pulse 
also strongly depends on the resulting Lorentz factor F t j. Thus it is possible to pro- 
duce a larger dynamic range of spectral peak energies by increasing T max , while 
leaving T m i n unchanged (see banded region in Figure 5). We found that more vari- 
able profiles have higher spectral peaks and can reasonably fit the observational 
data (Lloyd-Ronning & Ramirez-Ruiz 2002) provided that ( > 10" 2 . In this case, 
the synchrotron emission lies mainly below the BATSE window and the inverse 
Compton component carries most of the energy of burst. When synchrotron emis- 
sion dominates over inverse Compton scattering, as the latter takes place in the 
Klein Nishina regime (i.e. ( < 10" 2 ), we find that simply varying T max of the shells 
does not produced the dynamical range of peak energies that is observed. This is 
mainly because the typical observed synchrotron frequency, which at R c gives 

^ K — rJi+z) — ' (5) 



is unlikely to increase by simply incrementing r as R c oc T 2 gives hu syn oc T" 2 for 
a constant L w , e B , e e and F'^. Assuming the correlation between the Lorentz fac- 
tor and 9j tV , the isotropic energy itself is related to the source expansion velocity, 
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L w oc E 4n oc r max , and the typical observed synchrotron frequency is thenF] oc T" 3 / 2 . 
This may be overcome by implicitly assuming that the electrons are continually ac- 
celerated in the shocked region (see Lloyd & Petrosian 2000) or by postulating that 
the equipartition energy acquired by the electrons and the magnetic field (e e and e^) 
are dependent on T. While these are strong constraints in the emission properties, 
it is worthy of notice that the correlation between variability and spectral energy 
could be significantly affected by selection effects and more GRBs with indepen- 
dent spectroscopic redshifts are need it in order to test its validity (Lloyd-Ronning 
& Ramirez-Ruiz 2002). This relationship -if true- can help shed light on the rele- 
vant emission mechanisms in unsteady outflows. We would also like to emphasize 
that the wind optical thickness is one of the factors to which the variability and 
spectral peak energy is most sensitive. Indeed, it is reassuring that the inclusion of 
this effect is of the utmost importance to reproduce the power density spectrum of 
gamma-ray bursts arising from internal shocks in unsteady outflows (Panaitescu et 
al. 1999, Spada et al. 2000). 

3.4 Lethargic outflows 

Until this point we have assumed that the changes in the energy per unit solid 
angle are caused solely by an increase in the bulk Lorentz factor (M=constant). 
Nonetheless, it could be possible that changes in E 4w may be mainly caused by 
variations in the baryon loading of the wind. We have run such lethargic models and 
find that simply varying the baryon loading of the wind of the ejected shells does not 
reproduce the dynamic range of variability and spectral energies that is observed. 
We also study intermediate behaviors by assuming that the shell baryon loading 
and its Lorentz factor are related: T oc M a . We conclude that, within such models, 
the burst phenomenology can be accommodated provided that the changes in the 
energy per unit solid angle within different beams is dominated by T with a < 0.2. 
This is not surprising since the timescales that determine the pulse structure and 
thus the variability are highly dependent on R c oc Tj. An interesting consequence is 
that large variations in the energy per unit solid angle within different beam models 
should be strongly dominated by changes in V. 



4 Discussion 

The main conclusion we can draw from the simulations presented above is that, 
in the framework of the internal shocks model, there are realistic assumptions in 

4 This scaling is valid if eg, e e and are constant. However, it is possible to produce 
a larger dynamical range of spectral energies by increasing T max while leaving r m j n un- 
changed. In this case, T- • oc T and thus /iz/ syn oc T 1 / 2 
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which the variable activity of GRB is found to be correlated with absolute peak 
luminosity. We attribute this correlation to a variation of the energy per unit solid 
angle caused mainly by an increase in the emitting regions' Lorentz factors in the 
context of impulsive relativistic outflows. This change of the Doppler factor could 
be attributable either to a standard type of event viewed from different orientations 
or by changes in the collimation of GRBs jets. Within such models, collisions be- 
tween shells with higher Lorentz factors, which are intrinsically more luminous, 
produce more variable profiles. The existence of a correlation between the charac- 
teristic photon energy in the cosmological rest frame and the GRB variability, on 
the other hand, imposes more constraining requirements. The main characteristics 
of the modeled bursts with the above mentioned features are: a high electron injec- 
tion fraction - which implies that the inverse Compton scattering dominates over 
synchrotron emission - required to increase the spectral peak energy of the larger 
collision radii; and a wind optical thickness to scattering on cold electrons above 
unity is required to increase the duration of the pulses and decrease the photon en- 
ergy of the collisions that take place at smaller radii. 

We must also remain aware of other possibilities. For instance, we may be wrong 
in supposing that the main radiation mechanism assumed to be responsible for the 
burst event is synchrotron or its inverse Compton component. It could be that the 
GRB emission originates from a fireball moving out through an extended stellar 
envelope, along a funnel that is empty of matter but pervaded by thermal radiation 
from the funnel walls, leading to Compton drag (Lazzati et al. 1999; Ghisellini et al. 
2001). The fireball itself remains optically thick until it expands beyond the stellar 
surface (see also Ramirez-Ruiz, MacFadyen & Lazzati 2002). A burst with complex 
time- structure could then be produced by a series of expanding shells. In this sce- 
nario, faster moving shells will produce greater peak energies, as the burst emission 
simply reflects the temperature of the funnel photons up- scattered by the square of 
the bulk Lorentz factor. Furthermore, these fast shells will produce both smaller and 
more luminous pulses. This is because the observed variability timescale, which is 
related to the typical size of the region containing the dense seed photon divided 
by the time compression factor, increases with T and the Compton drag is more ef- 
ficient for faster shells. To this end, ejection of very rapid shells will lead to bursts 
produced by the ejection of very rapid shells are both more variable and have larger 
peak spectral energies. 

Much progress has been made in understanding how afterglows can arise from a 
forward shock or blast-wave moving into the external medium ahead of the ejecta, 
and in deriving the generic properties of the long wavelength afterglows that follow 
from this (Meszaros 2001). There still remain a number of mysteries, especially 
regarding the prompt emission, in particular, the formation of the ultra relativistic 
outflow, its structure and the radiation process. As we have shown, the correlations 
between the characteristic photon energy, the gamma-ray burst variability and its 
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luminosity, as well as the jet collimation, can help shed light on all three of these 
issues. 
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a) b) 




Fig. 1 . Diagram illustrating two different configurations of beaming models. The GRB phe- 
nomenology could be attributable either to a standard type of event with an axisymmetric 
energy distribution of the form E(0 V ) oc 6~P viewed from different orientations (a) or to a 
standard energy reservoir but with a variety of beaming angles (b). 
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Fig. 2. Evolution of the pulse duration as a function of the radius of collision and of the 
timescales that contribute to it. For lower Lorentz factor T max rs 10 2 the collisions take 
place at smaller radius where the wind is optically thick, photons are down scattered by the 
cold electrons before they escape, leading to an increase of the pulse duration A7d ~ A 7b 
and a decrease of the pulse energy. For T max ss 10 3 must collisions occur where the shells 
are optically thin A7d <C A 7b. The radiative cooling time is negligible with respect to AT a 
and ATa for collisions occurring at small radii, while for very large radii (i.e. small number 
of collisions) AT r is the dominant contribution. 
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Fig. 3. The luminosity per unit solid angle as a function of variability for bursts arising 
from multiple shock in a relativistic wind. The filled circles are bursts with secure redshifts 
estimates (Fenimore & Ramirez-Ruiz 2002a). The shaded area and banded areas represent 
the la regions of two different sets of random simulations.The former shows the effect of 
increasing r max /r m j n , while the latter depicts the effect of increasing T max but leaving r m ; n 
unchanged. The observed trend reproduced as bursts with higher source expansion velocity, 
which are intrinsically more luminous, produce more variable temporal profiles. 
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Fig. 4. Downscattered synchrotron and inverse Compton energy peaks as a function of the 
radius of collision. The dash lines represent the evolution of the spectral peak of these two 
emission components before downscattering (e e = 0.25, eg = 0.1, and C = !)• For collisions 
taking place at smaller radii (left panel), the wind is optically thick and photons are down- 
scattered by the cold electrons before they escape the wind. 
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Fig. 5. The peak of the vF v spectrum as a function of variability for bursts arising from an 
unsteady relativistic wind. The filled circles are bursts with secure redshifts estimates, while 
the empty circles are bursts in which the redshift is derived using the variability-luminosity 
distance indicator (Lloyd-Ronning & Ramirez-Ruiz 2002). The shaded and banded areas 
represent the la regions corresponding to the internal shock parameters of the simulations 
shown in Figure 3. Again, the observed trend is reproduced, as bursts with higher source ex- 
pansion velocity, which are intrinsically more luminous, radiate photons with higher char- 
acteristic energies after diffusion through the optical wind has taken place. 
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